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FRACTAL DIMENSIONS OF ROUGH DIFFERENTIAL EQUATIONS DRIVEN BY 
FRACTIONAL BROWNIAN MOTIONS 

SHUWEN LOU AND CHENG OUYANG 


ABSTRACT. In this work we study fractal properties of rough differential equations driven by a 
fractional Brownian motions with Hurst parameter H > i. In particular, we show that the Hausdorff 
dimension of the sample paths of the solution is min{d, -A } and that the Hausdorff dimension of 
the level set L x = {t G [e, 1] : X t = x} is 1 — dH with positive probability when d < -A. 
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1. Introduction 

Random dynamical systems are well established modeling tools for a variety of natural phenom¬ 
ena ranging from physics (fundamental and phenomenological) to chemistry and more recently to 
biology, economics, engineering sciences and mathematical finance. In many interesting models 
the lack of any regularity of the external inputs of the differential equation as functions of time is a 
technical difficulty that hampers their mathematical analysis. The theory of rough paths has been 
initially developed by T. Lyons lil2l in the 1990's to provide a framework to analyze a large class 
of driven differential equations and the precise relations between the driving signal and the output 
(that is the state, as function of time, of the controlled system). 

Rough paths theory provides a nice framework to study differential equations driven by Gaussian 
processes (see @). In particular, using rough paths theory, we may define solutions of stochastic 
differential equations driven by a fractional Brownian motion. Let us then consider 

rt d nt 

(1-1) X t = x+ / Vo(X s )ds + V / Vi(X s )dBi, 

Jo i=1 4o 
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where x € M n , Vq,V\,- ■■ , Vj are bounded smooth vector fields on and (B t )t >o is a d- 
dimensional fractional Brownian motion with Hurst parameter // G ( j. 1). Existence and unique¬ 
ness of solutions to the above equation can be found, for example, in |[T3l . In particular, when 
H = 4, this notion of solution coincides with the solution of the corresponding Stratonovitch sto¬ 
chastic differential equation. It is also clear now (cf. 0], [3] r4] [9J ) that under Hormander’s condition 
the law of the solution X t has a smooth density pt(x, y) with respect to the Lebesgue measure on 
R n . 

In this paper, we study fractal properties of the sample paths of solution X to equation (11.11) . 
More specifically, we investigate the Hausdorff dimensions of sample paths Af([0,1]) and level sets 
L x = {t € [0, l\,X t = x} of X on [0,1]. Our main result extends the classical result for fractional 
Brownian motions (see e.g. iflfill and ifTTl ) and is summarized as follows. 

Assumption 1.1. There exists a strictly positive constant A such that 

v*V(x)V*(x)v > A|n| 2 , for all v,x € M n , 

where we have set V = {Vj)i=i...., n ;j=i,...d- In particular, under this assumption, n = d in 
equation £3 and V 7 ],.... Vi are vector fields on M d . 

Theorem 1.2. Let X be the solution to equation ( 17.71 ). We have almost surely 


Moreover, for any given x € R d , 

i. If dH > 1, L x = 0 a.s.; 

ii. If dH < 1, dim% L x = 1 — dII with positive probability. 

Let us mention that some relevant results regarding sample paths properties of X has been 
studied in |2| under the same uniform ellipticity condition. In the aforementioned paper the authors 
have shown that for any bounded Borel set AcM^ 

P{X 4 hits A for t € [a, b ]} <(=^ Cap a (A) > 0 , 

where Cap a {A) is the a-dimensional Newtonian capacity of A. In particular, this characterization 
of hitting probability in terms of capacity implies that jj is the critical dimension concerning 
whether or not the process X hits a given point x in M d ; that is, X does not hit x almost surely if 
d > and X hits x with positive probability if d < -jj. The main ingredient in obtaining these 
results is upper and lower bound estimates for densities of X. 

Based on techniques developed in O, we are able to have some further study for density 
functions of X. In particular', we slightly improve the density estimate for the random vector 
(X s ,X t — X s ) to an exponential decay (see Theorem l3.4l belowl 

,, c f |6I 27 1 

< (t _ s y lH exp | C \ t _ s|272 j , 

for 7 < H. Let us point out that we will prove in Theorem 13. 3 1 help w that the density p/-s(f) of 
X t — X s has a faster decay rate in £, 

C f |£|(2tf+l)A2 A 


dim% AT([0,1]) = min < d , 
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This is due to the fact that we need to perform integration by parts twice in order to obtain the 
density of (X s ,Xt — X s ). For this purpose we have to sacrifice the decay rate in order to have an 
extra order of smoothness in the Malliavin sense. 

The organization of the paper is as follows. In Section 2, we give necessary preliminaries on 
rough paths and Malliavin calculus. In Section 3, we show several tail and density estimates for X 
that will be needed later in order to obtain our main result. Section 4 is then devoted to the proof 
of our main result. 


2. PRELIMINARY MATEIRAL 

For some fixed H > \, we consider (Q, T . P) the canonical probability space associated with 
the fractional Brownian motion (in short fBm) with Hurst parameter H. That is, Q = Co([0,1]) is 
the Banach space of continuous functions vanishing at zero equipped with the supremum norm, T 
is the Borel sigma-algebra and R is the unique probability measure on Q such that the canonical 
process B = {B t = (Bf ,..., Bf), t € [0,1]} is a fractional Brownian motion with Hurst 
parameter H. In this context, let us recall that B is a d-dimensional centered Gaussian process, 
whose covariance structure is induced by 

(2.1) R (t. s ) := E B J S Bf = i (s 2H + t 2H — \t — s \ 2H ) , s,t € [0, 1] and j = 1, ..., d. 

In particular it can be shown, by a standard application of Kolmogorov’s criterion, that B admits a 
continuous version whose paths are 7 -Holder continuous for any 7 < //. 


2.1. Rough paths above fractional Brownian motions. In this section, we recall some basic 
results in rough paths theory and how a fractional Brownian motion is lifted to be a rough path. 
More details can be found in Q and fl3l . For N € N, recall that the truncated algebra T N (M. d ) is 
defined by 

N 

T N (R d ) = 

771=0 

with the convention (R d )® 0 = R. The set T N ( R d ) is equipped with a straightforward vector 
space structure plus an multiplication <8- Let 7 r m be the projection on the rn-th tensor level. Then 
(T N (R d ), +, (8) is an associative algebra with unit element 1 € (R rf )®°. 

For s < t and m > 2, consider the simplex A™ = {(ui,..., u m ) € [s, t ] m ; u\ < ■ ■ ■ < u m }, 
while the simplices over [0, 1] will be denoted by A m . A continuous map x : A 2 —>• T A (W 1 ) 
is called a multiplicative functional if for s < u < t one has x s i = x SjU (8 x„./,. An important 
example arises from considering paths x with finite variation: for 0 < s < t we set 


x 


m _ 

s,t 




eq <8 ■ ■ • <8 e irn , 


where {e\,..., e,i } denotes the canonical basis of R rf , and then define the truncated signature of x 
as 

N 

Sn{x) : A 2 -A T n ( R d ), (s,t) S N (x) S:t := 1 + ^ x™ t . 

m=l 

The function Sn(x) for a smooth function x will be our typical example of multiplicative func¬ 
tional. Let us stress the fact that those elements take values in the strict subset G'- V (R' / ) C T N (R d ), 
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called free nilpotent group of step N, and is equipped with the classical Carnot-Caratheodory norm 
which we simply denote by | • |. For a path x G C([0, 1], G N the p-variation norm of x is 
defined to be 

ll x llp-var;[0,l] = Slip V |:Xj" 1 <g> X ti+1 | P 

nc[0,i] V i J 

where the supremum is taken over all subdivisions II of [0,1], 

With these notions in hand, let us briefly define what we mean by geometric rough path (we 
refer to (7} DEI for a complete overview): for p > 1, an element x : [ 0 , 1 ] —>• G I p (R d ) is said to be 
a geometric rough path if it is the p-var limit of a sequence (x rn ). In particular, it is an element 
of the space 

C p— var;[04] ([0)1]jG W (M d )) = {xGC([ 0 jl]jG W (M d )) : ||x|| p _ var;[0jl] < (X)}. 

According to the considerations above, in order to prove that a lift of a d-dimensional fBm as 
a geometric rough path exists it is sufficient to build enough iterated integrals of B by a limiting 
procedure. Towards this aim, a lot of the information concerning B is encoded in the rectangular 
increments of the covariance function R (defined by (12. 1 I f), which are given by 

Rt E [(B l t - B]) (.Bl - Bl)\ . 

We then call 2-dimensional p -variation of It the quantity 



/ 

( 

\ i Ip 

V P (R) = sup < 

1 \ ' +1 

1 \ Us i S i+l 

P ) ; ( s i)> (tj) e n 


\ 

< 

J 


where II stands again for the set of partitions of [0,1]. It is know that (see, for example Q) 
if a process has a covariance function with finite p-variation for p € [1, 2), it admits a lift to a 
geometric p-rough path for all p > 2p. As a consequence, we have the following for fractional 
Brownian motions: 

Proposition 2.1. For a fractional Brownian motion with Hurst parameter H, we have V p (R) < oo 
for all p > fjj. Consequently, for H > | the process B admits a lift B as a geometric rough path 
of order pfor any p > jj. 

2.2. Malliavin calculus for fractional Brownian motions. We introduce the basic framework of 
Malliavin calculus in this subsection. The reader is invited to read the corresponding chapters in 
lfl4l for further details. Let £ be the space of M . d -valued step functions on [0,1], and H the closure 
of £ for the scalar product: 

d 

((l[0,ti]J ■ ' ' ) l[0,td])) (l[0,si]> ; Rjtjj s i)- 

i= 1 

Let ei,..., ed be the canonical basis of W 1 , there is an isometry Kf : FI —> L 2 (\ 0,1]) such that 

k h( l[o,t] e i ) = ![o ,t] K H{t, •) ei, 
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where the kernel A'// is given by 


1 ft 3 1 

Kji(t,s) = CHS 2 ~ h (u - s) H -zu H ~i du, if > 1/2 

J S 

K H (t, S ) = CH ,1 (|) ^ (t- s) H ~ l > 2 + CH ,2 S 1/2 ~ H j (u- s) H ~^u H ~i du, H < 1/2 

J S 

for some constants ch, c rr. i • and ch, 2- 

Let us remark that TL is the reproducing kernel Hilbert space for B. Denote by M 1 the Cameron- 
Martin space of B, one can show that the operator 1Z := TZh '■ TL —>• M' given by 


U^-.= [ K H (; S )[K* H TP\(s)ds 
Jo 

defines an isometry between TL and Jtf. 

A A-measurable real valued random variable F is said to be cylindrical if it can be written, for 
a given n > 1, as 

F = f (B^ 1 ), .. ., B(r)) = f(j WldBs),..., dB s >) , 

where </>* € TL and / : K" —> Iff is a C°° bounded function with bounded derivatives. The set of 
cylindrical random variables is denoted S. 

The Malliavin derivative is defined as follows: for F € S, the derivative of F is the lfl' / valued 
stochastic process (D f F)o<t<i given by 


D t F = J2<P(t) 

i— 1 


K 

dxi 


(BW 1 ),..., B(<l> n )). 


More generally, we can introduce iterated derivatives. If F € <S, we set 


K,.,t k F = D fl 


D tk F. 


For any p > 1, it can be checked that the operator D k is closable from S into L p (fil-,TL® k ). We 
denote by D fc,p the closure of the class of cylindrical random variables with respect to the norm 




and 


E(F p ) + (ll DJF 

3 =1 

d°° = p| ^ d k ’ p . 

p> 1 k >1 



Definition 2.2. Let F = (F 1 ,..., F n ) be a random vector whose components are in ID) 00 . Define 
the Malliavin matrix of F by 

7f = ({DF,DF) w ) 

1 <i,j<n- 

Then F is called non-degenerate if "if is invertible a.s. and 

(det7 F ) -1 € n p >iL P (D). 
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It is a classical result that the law of a non-degenerate random vector F = (F 1 ...., F n ) admits 
a smooth density with respect to the Lebesgue measure on M n . Furthermore, denote by C^°(W l ) 
the space of smooth functions whose derivatives together with itself have polynomial growth. The 
following integration by parts formula allows to get more quantitative estimates: 

Proposition 2.3. LetF = (F 1 , F n ) be a non-degenerate random vector whose components are 
in B°°, and 7 p the Malliavin matrix ofF. Let G £ D°° and p be a function in the space Cjj°(M n ). 
Then for any multi-index a £ {1, 2,..., n} k , k > 1, there exists an element H a £ D°° such that 

E[dMF)G\ =E[<p(F)H a \. 


Moreover, the elements H a are recursively given by 

d 

%) = E^( G '(^)* JE) F J ) 

3 =1 

Hot — H(a k )(H( a i,...,a k _i))> 

and for l<p<q<oowe have 

whe re 1 = 1 + 1. 

Consider process W = { Wt , t € [0,1]} defined by 

W t = B((K* H )~\ l [0it] )). 

One can show that W is a Wiener process, and the process B has the integral representation 

B t = f K H {fs)dW s . 

J 0 

Based on the above representation, one can consider fractional Brownian motions and hence 
functionals of fractional Brownian motions as functionals of the underline Wiener process W. 
This observation allows one to perform Malliavin calculus with respect to the Wiener process W. 
We shall perform Malliavin calculus with respect to both B and W. In order to distinguish them, 
the Malliavin derivatives (and corresponding Sobolev spaces respectively) with respect to W will 
be denoted by D (and by D k,p respectively). The relation between the two operators D and D is 
given by the following (see e.g mi 

Proposition 2.4. Let D 1 2 be the Malliavin-Sobolev space corresponding to the Wiener process 
W. Then B 1,2 = D 1,2 and for any F £ F> 1,2 we have DF = FT*DF whenever both 

members of the relation are well defined. 


In order to estimate the bivariate density function for the random vector (X s ,X t ), we will need 
a version of conditional integration by parts formula. For this putpose, we choose to work on the 
underlying Wiener process W. The advantage of doing so is that projections on subspaces are 
easier to describe in a L 2 type setting. Set L 2 = L 2 ([t, 1]) and E/ = E(-|F/). For a random 
variable F and t £ [0,1], define for m > 0, p > 0, 
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and 

0.2) r Fli =(< D F<, D F.^) isi . £ii . 

The following formula is borrowed from lfl4l Proposition 2.1.4]: 


Proposition 2.5. Fix k > 1. Let F, Z s , G G (D°°) n be three random vectors where Z s G T s - 
measurable and (detr F+Zs ) _1 has finite moments of all orders. Let g G C£°(M d ). Then, for any 
multi-index a = (an, ..., a*,) G {1, ..., n} k , there exists a r.v. H^{F, G ) G n p >i n m >o D m,p 
such that 

IE. [(d a g)(F + Z s ) G] = E s \g(F + Z s ) H s a (F, G)|] , 
where H ( F. G) is recursively defined by 



H s a (F,G) = H s {ak) (F,H { 




(F,G)). 


Here 5 S denotes the Skorohod integral with respect to the Wiener process W on the inten’al [s, 1]. 
Furthermore, the following norm estimates hold true: 


\H s a (F,G)\\ p , s <c p jT^ s DF 


| k 

' k,2 k ~ 1 r,s 


\\G\\ k k , q , s , 


where - = - + -. 


3. Tail and density estimates 


Consider the following stochastic differential equation driven by a fractional Brownian motion 
with Hurst parameter H > 

(3.1) X t = x + V f VfXfidBl + f V 0 (X s )ds. 

i=i 4o 

Here V t , i = 0. 1,.... d are C°° -bounded vector fiefs on IP :" 1 which form a uniform elliptic system 
(recall Assumption ll.il ). Proposition 12 . ll ensures the existence of a lift of B as a geometrical rough 
path, which allows one to consider equation ( 13 . 11 ) as a rough differential equation. Then general 
rough paths theory (see e.g. J7] [8]]) together with some integrability results ((4] [51) shows that 
equation ( 13 . 11 ) admits a unique finite p-var continuous solution X in the rough paths sense, for any 
p > jj. Moreover, there exists A > 0 such that 


E 


exp A 


sup |^|^ +1 ) A2 
te[ o,t] 


< oo. 


The above moment estimate (or rather tail estimate) extends to the increments of X. 


Proposition 3.1. Assume the uniform elliptic condition for V \..... V,[. There exist non negative 
constants C and c not depending on s,t G [0,1] and £ such that 


sup 

s<u<v<t 


\x v - X u | > £ ) < Cexp 


c ^(2L7+l)A2 
(t — s) 2H 
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Proof. We borrow the idea from 0. Consider first the case 1/4 < H < 1/2. Taking up the 
notation of |0] we consider p > 2p and the control 

Wr n(u, V) = ||B|| P r 

’ / II Up—var; [u,v\ 

Then 0 Lemma 10.7] states that there exists constant cy only depends on the vector fields V's 
such that 

(3-2) \\X || p _var;[u,D] < C V (l|B|| p _ var; [ U)t ,] V ||B|| p _ var .j u ^ 

= Cy (u B ,p(u, v) 1/p V UJB,p(u, v)J . 

In particular, for any ti < t l+ \ we have 

(3-3) \SX ti t i+1 \ < cy ^[cUB,p(tj,tj+l)] 1/p V tj+l)) • 

Consider now a > 1 and construct a partition of [s, t] inductively in the following way: we set 
to = s and 

1 = inf < u > tp. ||B|| P r , i > a > . 

II Up—var;[q,«| — J 

We then set N aiS j tP = min{n > 0: t n > t}. Observe that since we have taken a > 1, inequal¬ 
ity (13.31) becomes 

\SX uti+1 \ < cy UJB,p(ti,ti+i) = cy a. 

Now for any fixed u < v in [s, t], suppose l and m are such that u € [tpti + 1 ] and v € [t m ,t m + 1 ]. 
We have 


\x» -Xv\<Yl \ dx ut i+1 \ + \^ tl ,u\ + \SX uM+1 \ 

i =0 

TTL Na,s,t,p 

+ E !“«.«! + l«wl + + E \6X titi+1 \ 

i=l- 1-1 i=m +1 

<Cy a (N at t,p + 2 ) 

Since fractional Brownian motions have stationary increments, Theorem 6.4 in JH implies that 
there exists constants C and c not depending on s, t and £ such that 

r ti H+l 

( 3 . 4 ) i 4 

This easily yields 


f ct 2H+l 1 

P (X a ,s,t,p + 2 > 0 < C exp | ~ ^ _ s y H j 


( < SUP <t ^ ~~ Xu][ - ^) - V ( Cva ( N a,t,P + 2) > 0 < C exp | — 2if |, 


which is our claim. The case H > 1/2 is handled along the same lines, except that the coefficient 
£ 2 H+i j n (|3 4 |) j s replaced by £ 2 . □ 


The vector X t is a typical example of a smooth random variable in the Malliavin sense. Recall 
that D is the Malliavin derivative operator with respect to the underline Wiener process W and T p,t 
is defined in (12.21 ) for a random variable F. The following estimate is a restatement of Proposition 
5.9 in El. 
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Lemma 3.2. Let e € (0,1), and consider H € (1/4,1). Then there exist constants C,r > 0 
depending on e such that for e < s < t < 1 the following holds: 


lip —1 lid 

II 1 X t -X s ,s\\d,2 d + 2 ,s 


\\D(X t 


X s)\\d,2 d + 2 ,s 


< 


c 


0 t-s ) 


2 dH 


Ef +2 (1 + G) 


< C(t-s) 


dH E? d+2 


(1 + G), 


where G is a random variable smooth in the Malliavin sense and has finite moments to any order. 


With the above le mm a in hand, we are able to obtain an upper bound for the density of X t — X s . 


Theorem 3.3. Fix e € (0,1). Let pt— s (z) be the density function of Xt — X s for s,t € [e, 1]. There 
exits a positive constant C depending on e such that for all s < t, we have 

C ( M( 2H+1)A2 1 

Proof We first write 


p t - s {z) = E[d z (X t -X s )\ 

= EE S [8 z (X t - X,)]. 


Next, we bound M st = E s [5 z (X t — X s )] by first using the conditional integration by parts formula 
in Proposition 12. 5 1 and then Cauchy-Schwarz inequality. 


< C\\T x ]_ x J\ d d ^ !S || D(X t - X s )\\l 2d+2 s E V 2 [l (Xt -x s > 6 )] • 

Thus, owing to Lemma [T2l we obtain: 

C 


(3-5) 


Pt-s(z) < 


(t - s) 


dH 


E 


E 


d 

2 d + 1 


(1 + G) 2 E \ /2 [l(x*-x s >a)] 


Then an application of Cauchy-Schwarz inequality together with Proposition ^, ll finishes the proof. 

□ 


Finally we slightly improve the estimate in |j2]| for the bivariate density of (X s , Xf). 

Theorem 3.4. Fix e € (0, 1) and 7 < H. Let p s ,t(zi, z 'i) be the joint density of the random vector 
(X s , Xt) for s,t € [e, 1]. There is a positive constant C depending on e such that for all s < t, we 
have 

t ■, / C J ki- z 2 | 27 1 

Ps,t( z uz 2 ) < _ g)dH exp | _ s|272 I • 

Proof. The proof is similar to the one in the above, but with a slightly more careful estimate for 
the tail probability of X t — X s . 

First note the existence and smoothness of the density function p s> t(zi,Z 2 ) a consequence of 
Proposition 13.21 We then write 

Pa,t(zi,z 2 ) = Ps,t-s{zi,Z 2 - 21 ), for Zl,Z 2 € M'/ 










10 


SHUWEN LOU AND CHENG OUYANG 


where p s>t _ s (-, •) denotes the density of the random vector (X s , X t — X s ). We now bound the 
function p s ,t-s, which shall be expressed as 


= E [% (X s ) 5 ? 2 (X t -x s )}, for 6 , 6 € R d , 

= E[%(X S )E S [d^ 2 (X t - X s )}} . 

As before, we can bound M st = E s [d^ 2 (X t — X s )] as follows. 

\M s ,t\ < C\\V- x ]_ x J\ d dfid+2jS || D(X t - X s )\\ d d ^A ' 2 [ 1 {x t -x.>b)\ 

Thus, owing to Lemma [3721 we obtain: 

(3.6) p s ,t-s(Z 1,6) < {t _ C s) nH E %(^)E|^ rT (l + G) 2 EV 2 [l {Xt _ Xs>€2) ] 

To proceed, recall that B is the lift of B as a rough path. Set 

H rl d(B v ,B u )<i 


A^y, g (B) = 



/0 Jo \v-U |T"? 


dudv, 


where 7 < H and q > 0. The reason we need to consider these Besov norms As^(B) is that they 
are smooth in the Malliavin sense. From the Besov-Holder embedding we have for any e > 0 such 
that 7 + e < H, and large enough q 

||B||o,i ;7 < ca/Vk^B) 1 ^. 

Hence, it is readily checked by (13.21) that (note that 0<f — s<lin our case) 

\Xt-X s \<C\t- s| 7 (l + A/ 7 + ei g(B) 1 / 9 ) 1 / 7 . 

Furthermore, by [| 6 j], for any 7 < H there exists a constant A > 0 such that 

Eexp {A||B||o 1 7 } < 00 . 


This together with 

A/’ 7 + e ,q(B) 1//<? < ||B|| 0l l, 7 + e , 

implies that A 7 + f , f; (B ) l/ ' / has Gaussian tail. Thus, for large enough q and small enough A, we 
have for some constant C > 0 

Es [1(W-W>6)] < ^exp {-A } E s exp {a(1 + M 1+t , q (B))^} . 

Plugging this inequality into (13.61) . we end up with: 

( 3 . 7 ) P, f -,((ub) < 7Z7777 ex P 

where T'i and 'I 7 arc two random variables which are smooth in the Malliavin calculus sense. 
Also note that we may need to choose A even smaller to make sure Malliavin derivatives of <J >2 
have moments to certain large order. Based on the above consideration, we can now integrate (13.71) 
safely by parts in order to regularize the term % (X s ), which finishes the proof. 

□ 


To close the discussion in this section, let us state a mild lower bound (strict positivity) for the 
density of X t . We direct interested reader to (3 for its proof. 
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Theorem 3.5. Consider the solution X to equation (ED driven by a d-dimensional fractional 
Brownian motion with Hurst parameter H > j. For each fixed x G W l , denote by pt(x, •) : R d — >• 
M_l_ the density function of the random variable X/ . Assume that Assumption [7771 is satisfied. Then 
Pt(x , y) > Ofor all y G R d . 


4. Fractal dimensions of SDEs driven by fractional Brownian motions 


We first briefly recall the definition of capacity and packing dimension, as well as their connec¬ 
tion to Hausdorff dimension. A kernel k is a measurable function k : R d x R d —>■ [0, oo]. For a 
Borel measure p on R d , the energy of p with respect to the kernel n is defined by 





K(x,y)p(dx)p(dy). 


For any Borel set E C R d , the capacity of E with respect to n, denoted by C K {E), is defined by 

r -i-l 


c K (E) 


inf 

MG V{E) 




Flere V(E) is the family of probability measures carried by E. Note that C K (E) > 0 if and only 
if there is a probability measure p on E with finite K-cncrgy. Throughout our discussion, we will 
mostly consider the case when k(x, y) = /( \x — y |), where 


fir) 


r a , if a > 0; 

lo S (?Xi) > if a = 0. 


The corresponding energy and capacity will be denoted by £ a {p) and C a (E), respectively; and the 
former will be called the a-energy of p and the latter will be called the Bessel-Riesz capacity of E 
of order a. The capacity dimension of E is defined by 


dim c (£') = sup{a > 0 ]C a {E) > 0}. 

It is know by Frostman’s theorem (cf. fflOl or HD) that 


dim% E = dim c (£ 1 ), 

for every compact subset E of W 1 . Hence, in order to show dim-^ E > a one only needs to find a 
measure p on E such that the a-energy of p is finite. 

Packing dimension and packing measure were introduced as dual concept to Hausdorff dimen¬ 
sion and Hausdorff measure. Denote by dimp by packing dimension. For any e > 0 and any 
bounded set F c R d , let N(F, e) be the smallest number of balls of radius e (in Euclidean metric) 
needed to cover F. Then the upper box-counting dimension of F is 


dim is F 


log N(F,e) 
lim sup---. 

ej .0 -loge 


The packing dimension of F can be defined by 


dim-p F = inf < supdim^F n : / C U™ =1 F n 

{ n 


It is known that for any bounded set F C R d (cf. lUTTO . 


(4.1) 


dim-p F < dimp F < dimgF 1 < d. 
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4.1. Hausdorff dimension for SDE driven by fBm. Recall that X is the solution to 

rt d rt 

(4.2) X t = x+ / Vo(X s )ds + V / Vi(X s )dBi, 

Jo i=1 Jo 

where V, i = 0,1,..., d are C°° -bounded vector fiefs on W 1 satisfying Assumption ll.il Denote by 

X([0,l]) = {X(t):te[0,l}} 
the collection of sample paths of X on time interval [0,1]. 

Theorem 4.1. We have 

dim B X([ 0 , 1 ]) < a.s. 


Proof. For any constant 0 < 7 < H, we have 

\X t - X s 


(4.3) 


sup 

s,ie[o,i] 


|s - t\l 


= \\X\ 


7; [0,1]) 


and \\X || 7 ; [o,i] has moments to any order. We fix a sample path w and suppress it. For any integer 
n > 2, we divide [0,1] into m n sub-intervals {R n ^} with length n~ l ^ H . Then 


m n < ciuh 


and X([0,1]) can be covered by X{R Ui i ), 1 < i < m n . By (14.31) . we see that the diameter of the 
image X(R. n i) is controlled by 

diamX (Rn,i) < C 2 n~ 1+S , 

where 5 = (H — 7 )/H. Consequently, for e n = C 2 n” 1+<5 , X([0,1]) can be covered by at most m n 
balls in of radius e n . That is 

N(X([0 ,1 },£ n ) < cinTi. 


This implies 

dim B W([ 0 , 1 ]) < a - s - 

Letting 7 f H finishes the proof. 


□ 


Next, we turn to the lower bound. 

Theorem 4.2. Let X be the solution to equation ( IJ.7D . We have 

Proof. Note that for any e > 0, we have dim% X([0.1]) > dim% X(\e, 1]). To prove the claimed 
result, it suffices to show that 

dim-H X([e, 1]) > 7 a.s. 


dim^ X ([ 0 , 1 ]) > min < d 


for every 0 < 7 < rriin{ri. 1/H}. 
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Let f ix be the image of measure of the Lebesgue measure on [e. 1] under the mapping t i->- X t . 
Then the energy of px of order 7 can be written as 

px(dx)px(dy) f f dsdt 


/ / 

it" JW 


\x-y\~ 


L 


[e,l] J[e, 1] I X t - X sV 


Hence, by Frostman’s theorem, it is sufficient to show that for every 0 < 7 < min{d, 1 /H}, 


(4.4) 

By Theorem 13.31 



|*|(2H+1)A2 
e C\t-s\ 2H d z 

C f 1 \ |^|(2H+l)A2 

“ \t-s\^ H J Rd jffr 6XP |£ — s |C[2i/-(2i/+l)A2] 

a 

- |f - sY> H ' 


Et 


\Xt-X a \r 


< 


/ n 

J R d \ z \ 


c 


Ft-s 


I dH 


di 


In the above, we have used the fact that 7 < d to ensure the integration near 0 is finite. We have 
also used that 2 H — (2II + 1) A 2 > Oto ensure the integration at infinity is finite and uniform 
in s, t € [0,1]. Now plug the above into (14.41) and note that we assume 7 < -jj. The proof is thus 
completed. □ 


Theorem 14.11 and Theorem 14.21 together with (14.11) give us the Hausdorff dimension of sample 
path of X. 


Theorem 4.3. Let X be the solution to equation ( 14.21 ). We have almost surely 


dim?* X( [0,1]) 


min 



It is expected that similar argument also allows us to analyze the Hausdorff dimension of 
GrX([0,l]) = {(t,X(t)) : t € [0,1]}, the graph of X on [0,1]. The result is summarized in 
the following theorem. 


Theorem 4.4. The following holds regarding the dimension of the graph of X. 


dim GrX([0,1]) = min | (1 — H)d + 1 , 


i.e., 


dimGr(X[ 0 , 1 ]) = 


Hi 


Hd > 1; 
l-H)d + l, Hd< 1 


Proof We first prove the upper bound of dim Gr(X[0,1]). Given any <) > 0, we can cover 
GrX([0,1]) by m n many balls in M 1+,i with radius n~ 1+s , where m n < cn l / H , which provides 
an upper bound l/H for dim Gr (X[0,1]). To show the other upper bound (1 — H)d + 1, we 











14 


SHUWEN LOU AND CHENG OUYANG 


notice that as in the proof to Theorem 14.11 each R n ^ x Xr u . can be covered by l n balls with radius 
n~ l / H , where 

l n = n (-5’+(- 1 +' 5 )) d+1 +i. 

Therefore Gr2f([0,1]) can be covered by m n x l n balls with radius n~ x l H . Letting 5 > 0 go to 
zero proves that 


dim GrX([0. 1 ]) < H 




1 + (1 - H)d. 


Now we need to prove the lower bound for dim Gr2f([0,1]). Since dim GrX > dim% X 
always holds, when (1 — H)d +1 > 1/H, or equivalently, Hd> 1, dim GrX > dim-^ X = l/H. 
Therefore we only need to worry about the case when (1 — H)d + 1 < 1 ///, or equivalently, 
Hd < 1. For 7 such that d < 1/H < 7 < 1 + (1 — H)d arbitrarily fixed, we will prove 
dim GrX([0,1]) a.s.. By Frostman’s theorem, we will show 


(4.5) 


f J 

Jt =0 Js =0 


E 


(\t-s\ 2 + \X t -X s \ 2 ) 


2\7/2 


dsdt < 00 . 


By the upper bound density estimate for X t — X s , we have for the left hand side of (14.51 ) that 


l-uL-n ' [<|t — s | 2 + \x t — x s \ 2 yi 2 

£ f f 

Ite[ 0 , 1 ] Jse[ 0 , 1 ] I y m d K — s\ dH+ ^ 


t =0 Js =0 JXt—X s :=ySM d 1^ — s \ dH 

1 


dsdt 


1 \y\( 2H + l)A2 

e 


(\t — s | 2 + \y\ 2 ) 


7/2 


dydsdt 


| J/ |(2H+l)A2 
- JB - 


< 


/ / 

Jte[ o,i] Js&[ o,i] 

/ / 

dtefoal dsi 


1 


If — s | 1-<5 


ie[0,l] JsG[0,l] 


f — s 


1-6 


L 

i 


e \t-s\ 2H dydsdt 

d -1 r (2H+l)A2 


r 

|f — s| d 

7 ,(2n+l)A2 

e \t-s\‘ 2H d 


e It-sF 3- drdsdt 


|f — s| 


dsdt < oo 


for any 5 > 0 such that 7 + dH < 1 — 6 + d. The integral converges since 2 H < (2 H + 1) A 2. 
Therefore it has been shown that when (1 — H)d+1 < 1/H, i.e, when Hd < 1, dimGrX([0,1]) > 
(1 — H)d + 1. This proves that 


dimGrX([ 0 , 1 ]) > min 


(i-H)d+i,±y 


□ 


4.2. Hausdorff dimension of level sets. In this section, we study Hausdorff dimension of the 
level sets of processes driven by fBM. Before we state the main results, we have the following two 
lemmas regarding the upper and lower bounds of Hausdorff dimensions of the level sets which 
actually matches each other. The following lemma addresses the upper bound. 
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Lemma 4.5. Let X = {X(t), t £ K} be a process driven byfBM on t e [e, 1]. 
(4.6) dimgLa; < 1 — dH , a.s. 

and L x = 0 when the right hand side of m is negative. 


Proof. For a fixed integer n > 1, we divide the interval [0,1] into m n subintervals m n 

with lengths n~ l ^ H . m n < cn l / H for some c > 0. Let 0 < 5 < 1 be fixed, and let ti be the 
left endpoint of In the following proof, c is always some strictly positive constant which may 
change from line to line. We have for 1 < i < m n 

(4.7) 

P(i€ X(Ii)) < P ( max \X 8 - X t \ < n“ (1_,5) ;x € Xilf) ) + P ( max \X S - X t \ > n~ {1_,5) 
\s,teh ) 



We first claim that the first term on the right hand side of (14.7b is bounded by cn *4 Assuming 
the process is started at z and applying Theorem 13.41 we have for 1 < i < m n , 


max |X S — Xf\ < n *4 s ^;x G X(If) 

s,teii 


< P (\X ti -x\< 

<-L 
<-L 


ly^B(x,n (- 1 <*)) ^ 


|y- g |(2g+l)A2 

dH e ctr dy 


y£B(x,n _ ( 1_lS )) £ 


]v _ x] (2H+l)A2 
C S 377 - J 

dw e 1 d y 


rn- C 1 -*) -dr-1 r^(2H + l) 

< c I —-—e dr < cn ~ lyl ~^ d , 


t r =0 


£ 


tiff 


where the last inequality is due to the fact that ti is bounded away from zero by at least e. 

For the second term on the right hand side of (14.7b . we would like to bound it from above by 
e~ cn , which will then bound the left hand side of (14.7b by cn~'' 1 ^ d . Indeed, it has been shown 


(4.8) 


It follows from (14.8b that for 1 < i < m n . 


sup \Xf — x\ > £ < cexp 

.se[o,t] / 


Ch ^2H +t)A2 


, for t € [ 0 , 1 ]. 


max \ X S — XA > n ^ ^ 


S,t£li 

Therefore, we have shown that 


„-(l-<5)[(2»+l)A2] 


< ce 


(2 H) 


< e 


,-2(1 S) 

- 


= e 


(x € X(Ii)) < P f max \X 8 - X t \ < n _ ( 1_a) ;x G X{If) J + P f max \X 8 - X t \ > n~( 1 ~ 5) ) 

J \^S )t£li J 

< cn~( l -^ d + e~ cn2S < cn~( 1 ~ s ' >d , 1 < i < m n 


(4.9) 


for sufficiently large n. 

If dH > 1, we choose <5 > 0 such that (1 — S)d > l/H. We denote by N n the number of 
times that X(/.,) visits x, as i increases from 1 to m n . It now follows from (14.9b and the fact that 
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m n < cn}^ H 

(4.10) E(JV n ) < cn 1/H • n~^- s)d -> 0, as n -> oo, 

which combined with Fatou’s lemma yields that A r n = 0 for infinitely many n, since the random 
variables N n are integer-valued. Therefore L x = 0 a.s.. 

| of L x can be defined by setting f = f if x € X(Ii), and U = 0 

otherwise. Each interval f has length n~ l ^ H . Let the number of such nonempty I, be denoted by 
M n . In the same way as getting (14.101) . again by (14.91 ) we have 

(4.11) E(Af n ) < cn- 1/H n-( 1 ~ s)d = cn 1/H “ (1_5)d . 

By picking 6 > 0 sufficiently small, we may define a constant rj by 

rj = 1/H — (1 — 26)d > 0. 

We consider the sequence of integers n 3 = 2 J . j > 1. (14.1 II) and Chebyshev’s inequality imply 

IP {M nj > e) < - 2 J (^ _(1 “ ,5)d ), Ve > 0 . 

Borell-Cantelli lemma yields almost surely M nj < cn - for all j large enough. This implies that 
dim gL x < Hrj almost surely. Letting 5 | 0 along rational numbers proves 

dim bL x < 1 — Hd a.s., 

which proves the upper bound. □ 

The next lemma is on the lower bound of Hausdorff dimension of L x . 


Lemma 4.6. There exists some C > 0 such that 

P{dim% L x > 1 — Hd} > C. 

Proof. To prove the lower bound for dim^ L x , we fix a small constant <5 > 0 such that 

7 := 1- (1 + 5)Hd > 0. 

Note that if we can prove that there is a constant c > 0 independent of 5 such that 
(4.12) P {dim w L x > 7 } > c, 

then the lower bound will follow by letting 5 | 0. The strategy to prove (14.121) is standard and 
based on the capacity argument by Kahane ifTOl . We spell out all the details for the convenience of 
the readers. 

Let be the space of all non-negative measures on M with finite 7 -energy. It is known that 
is a complete metric space under the metric || • || 7 given by 


fi(dt)ti(ds) 


\t - s|7 

We define a sequence of random positive measures //„. := r, ■) on the Borel sets C of [e, 1] by 

n\Xt — x\ 2 ' 

‘limy exp [ - 

>c 


(4.13) 


Tn{C) = f ( 2 vrn ) d/2 

Jc 


dt. 
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It follows from Kahane fflOil or Testard lfl5l that if there exist positive constants c,, ! < i < 3 such 
that 


(4.14) E(||/tt n ||) > ci, E (||/t n || 2 ) < c 2 , E(||//„|| 7 ) < c 3 for all n., 

where \\p n \\ = /i n ([e, 1]) denotes the total mass of fi n , then there is a subsequence of { u n }, say 
{ n nk }, such that /i„ fc —>• // 7 .M+, and // is strictly positive with probability at least cf/(2c2). It 
then follows from (14.131) that y has its support in L x almost surely. Moreover, the third inequality 
of (14.131 ) together with monotone convergence theorem imply that the 7— energy of y is finite. 
Therefore Frostman's theorem yields (14.12b with c = c 2 /(2o)). It remains to verify the three 
inequalities in (14.141) . 

In the following computation, the strictly positive constant c may change from line to line. 



In the last line above we have performed a change of variable y/n{y — x) = z — x and restricted 
the integration over z in a ball \z — x\ < 1. Then the desired low bound clearly follows by the fact 
that p(t, x, y) is globally strictly positive (Theorem [33]). 

For the second inequality of (14.141) . by applying the transition density estimate in Theorem 13.41 
we get 


E(||/tn|| 2 ) = E 


= 2 • (2nn) d ■ E 


't€[e,l] 


(2itn) d/2 


n I Xt — x\ 


exp 


dt 


[ (2nn) d/2 

J sE[e,l] 


n IX, — x\ 


exp 


ds 


f f 


exp 


n (\X t — x\ + \X S — x\ 


t=e J s=e 


<cn-f [ f 1 /"exp 

Jz&R d Jy£R d Jt=e Js=e 
rl ri 


11 i 2 1 1 1 2 

n [\y — x\ + \z — x\ 


dsdt 


<«-/ f j j 

J z£R d J y&R d Jt—e J s 


t=e J s=e 1 ^ s 


Hd 


exp 


|f — s\ Hd 

n (\y — x\ 2 + \z — x| 2 ^ 


dsdtdydz 


dsdtdydz 


= cn 


[ [ e- : * L r t dy = cn d - < c 

Jtd[e, 1] Jse[e,l] r s l j zd R d J yd R d 
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where the integration in t and s converges because Hd < 1 . 
Similarly, again by Theorem 13.41 


f f e 


e||a*X = 2 


7 t n) d exp 


n (\X t - x\ 2 + \X S - x \ 2 ) 


|i - S |7 




f 1 r f 1 [ ( , sH 

n fix — y 2 + la; — z\ 2 ) 

dt ds / (27rn)“ exp 


Jt=e Js=e J z£ J y& M d 

2 


|i — s|7 \t — s 


HcL 


exp < —A 


I y - ^l 27 

'C\t- s| 2 7 2 


dydz 


f f 


< C 


jdsdt , 


lt=e Js=c I* ~ s\^ Hd 

which converges since 7 = 1 — (1 + 5)Hd. This proves the third inequality of (14.141) . The proof 
to the upper bound is thus finished. □ 


The following theorem addressing the Hausdorff dimension of the level sets of processes driven 
by fBM is an immediate consequence of the combination of Lemma 1431 and Lemma l4~6l 

Theorem 4.7. Let X = {X(t),t € M} be a process driven by fBM on t € [e, 1]. 

(i) If dH > 1, then for every x € M. d , L x = 0 a.s. 

(ii) If dH < 1, then for every x € 

dim% L x = dinEp L x = 1 — dH 

holds with positive probability. 
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